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1. ABSTRACT 

This paper describes the operations to control the Moon 
impact of the 3 -axis stabilized spacecraft SMART- 1 in 
September 2006. 

SMART- 1 was launched on 27/09/2003. It was the first 
ESA mission to use an Electric Propulsion (EP) engine 
as the main motor to spiral out of the Earth gravity field 
and reach a scientific moon orbit [1]. During September 
2005 the last EP manoeuvres were performed using the 
remaining Xenon, in order to compensate for the 3 rd 
body perturbations of the Sun and Earth. These 
operations extended the mission for an additional year. 
Afterwards the EP performance became unpredictable 
and low, so that no meaningful operation for the moon 
impact could be done. To move the predicted impact 
point on the 16/8/2006 into visibility from Earth an 
alternative Delta-V strategy was designed. Due to their 
alignment, the attitude thrusters could not be used 
directly to generate the Delta-V, so this strategy was 
based on controlled angular momentum biasing. Firing 
along the velocity vector around apolune, the remaining 
Hydrazine left from the attitude control budget was 
used, to shift the impact to the required coordinates. 

2. INTRODUCTION 


The SMART- 1 AOCS allows the S/C to follow any 
attitude profile commanded by the ground. The attitude 
control is based on reaction wheels as actuators for all 



Fig. 1: SMART- 1 with attitude thruster pairs at the 
bottom. Illustration by AOES Medialab, ESA 2002. 


nominal modes and IN attitude (Hydrazine) thrusters 
for reaction wheel desaturation and for rate reduction in 
case of contingencies. These 2x4 attitude thrusters are 
mounted along the Z-axis away from the centre of mass, 
pointing in the X/Y-plane of the S/C. This means that 
Delta-V ( AV ) and Delta-L (A L , where L is the S/C 
total angular momentum) can occur simultaneously. 
Two autonomous star trackers, 5 attitude rate sensors 
(gyroscopes) and 3 coarse sun sensors are utilized as 
attitude sensors. 

The EP engine is mounted on a two axes articulation 
mechanism (EPMEC) attached to the S/C. While the 
engine is firing, a control loop of the AOCS aims to 
reduce the total S/C angular momentum by rotating the 
EP engine with respect to the S/C-body. The control 
loop reaches a steady state when the EPMEC- 
articulation is such that the thrust direction points to the 
S/C centre of mass. In this steady state, the thrust 
direction is about 4 degrees away from its mid-position, 
in which the thrust vector is aligned with the S/C Z-axis. 
Typical excursions from the steady state rotation angles 
are of the order of 1 degree. As a result of this 
momentum management with the EP thrust, which 
worked better than expected, the Hydrazine budget was 
not exhausted, so that approximately 6 kg of Hydrazine 
were left, and could be used for the required A V at the 
end of the mission. 

The solar arrays (SA) are mounted along the S/C 7-axis, 
and can be rotated infinitely. The Star Trackers (ST) are 
mounted on the +/- 7-faces with their boresights about 
40 degrees away from the +/- 7-axis. Therefore, the Sun 
will never be in the field of view (FOV) of a ST as long 
as the SA panels are perpendicular to the Sun direction 
(S). However, it is possible for a ST to be blinded by the 
Earth or the Moon. Simultaneous blinding of both STs 
must be avoided, since the attitude rate sensors do not 
permit the outage of both STs for more than one hour. 
This was especially delicate at times when the S/C was 
at low moon altitude. Apart from constraining the 
attitude, the STs also impose a constraint on the S/C 
rate: To guarantee proper tracking the angular S/C rate 
was qualified to 0.15 degrees/s [2]. Due to a limitation 
in the torque that can be exerted on the reaction wheels, 
the angular acceleration of the S/C was qualified to 
0.0005 degree/s 2 . 
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2.1 Initial situation 

The scientific community requested the impact point 
over the Moon’s terminator, with the best illumination 
conditions for impact observation. This request was 
driven by various scientific objectives. Namely 
spectrometric information constraining the mineralogy, 
the thermal and dynamic evolution of the impact flash 
and the ejecta morphology [3]. 

After about 2 years in operation, SMART- 1 reached the 
operational Moon orbit in early 2005. Six months later 
in September 2005 a further series of EP manoeuvres 
was performed to extend the mission. These 
manoeuvres were symmetric to the perilune, so the line 
of apsides is rotated to a point where the Earth and Sun 
perturbations initially raise the perilune again. 

To extend the operational period to a maximum the 
second extension took into account that the Xenon was 
exhausting. Therefore, the strategy for the last firing 
period of the EP had to consider that it could be stopped 
at any time. In the end, almost all of the second firing 
was performed and from this time on the engine was 
declared no longer reliable (and operation was stopped). 
This left SMART- 1 in a free drift orbit from where the 
S/C would eventually crash into the Moon, due to the 
increasing eccentricity. In fact the orbit propagation 
predicted (for this free drift orbit), that the security 
height of 300 km would be violated on the 8th June 
2006 and the S/C would crash on the 17th August 2006. 
However since the prediction period was almost one 
year, these results had an error in the impact time of 
more than a day due to the fact that irregularities of the 
Moon’s surface where not precisely known. 

The first part of the paper describes how the Delta-V 
was implemented using the four attitude thrusters, the 


pointing profiles [4] required and the S/C dynamics 
involved [5]. The second part shows how the 
manoeuvres were calculated and optimised to achieve 
the scientific requirements of the impact. 

3. IMPLEMENTATION OF THE DELTA-V 

3.1 Chosen strategy 

To implement the required AV , the following idea for 
the usage of the attitude thrusters was chosen: 
Command alternating target angular momentum along 
the symmetry axis of the solar arrays, i.e. the y-axis of 
the S/C, with synchronically alternating attitude around 
the same axis. The parameters sent to the S/C are the 
target levels of the wheels. Since the attitude thrusters 
are asymmetric about the S/C centre of mass, this will 
lead to a AV , which can be used to adjust the orbit 
[Fig. 2] and [Fig. 3]. 



Time (hours) 

Fig. 3: Commanded quaternion profile: Apart from 
rotating the body, adjustments for the pointing direction 
and the constraint of sun aspect angle (to 90 degrees) 
are made [4] 



Fig. 2 : Chosen strategy: To achieve the required Delta-V seven Delta-L manoeuvres per orbit are made with re- 
orientation of the S/C body and therefore the thrusters in-between [6] 
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This option was preferred amongst others, since the S/C 
had many constraints, the efficiency was high, the 
required performance could be reached and the 
dynamics was controllable. 

3.2 RCS model and AOCS controller 

Each Delta-L will linearly translate into a Delta-V : 

^ WoL ~ M force Torque ^ 

since A V WoL = M Forc \n anc * ^SL — M Torque T on , 

where M Force , M Torqu ~ X and 7;„are the force matrix, 

the pseudo inverse of the torque matrix and the thruster 
on-times respectively. 

Apart from compensation for the accumulated 
disturbance torques, the first A L will load 2 Nms along 
the Y-axis. Each subsequent A L (of 4 Nms) will flip 
this to -2 Nms along the same axis. Finally, a last A L 
will bring the total S/C angular momentum back to zero. 
Therefore, the resulting AV in the S/C frame will also 
alter in-between +/- 0.0264 m/s (or half of that) along 
the X-axis. The geometry of the three torque controlled 
reaction wheels, is such that one of then will only be 
used for the compensation of the attitude and external 
disturbances. 

Finally flipping the S/C body around the Y-axis by 180 
degrees after each A L , will cause the AV WoL in the 

inertial frame to accumulate in the required direction 
[Fig. 4]. These flips have two advantages: rotations take 
place around the smallest moment of inertia of the S/C 
and the rotations can be chosen such that the loaded 
angular momentum will temporarily decrease during 


each slew [Fig. 4], so that no extra margins or 
constraints occur. 

Using the techniques as described in [4] and applying an 
improved slew (with the difference that the angular 
velocity is now split in three phases: acceleration, coast 
arc and deceleration [Fig. 4]) these manoeuvres can be 
performed fulfilling all S/C operational constraints, 
namely keeping the SA panels perpendicular to the Sun. 

However, since the resulting A V W()L appeared too 

small, there were several options considered to improve 
this situation: 

• Usage of all four reaction wheels. This would 
allow to load twice as much A L . 

• Shorten the time for the slews by violating the 
rate constraint. 

• Speed up the duration of the A L by adjusting 
the duty cycle. 

After an investigation from ESOC/FD and SSC [5] 
option 2 was chosen, such that the maximal S/C rate 
was 0.23 degrees/second, since the limiting factors the 
controller of the SA mechanism and the Star Trackers 
proved to be able to cope with this rate. The angular 
acceleration of the S/C was set to 0.00375 degree/s 2 , 
well within the values of [5], but above initial 
qualification. This was possible due to the very solid 
AOCS controller of SMART- 1 and relatively low 
pointing accuracy requirements. It means that all 
standard operational constraints were fulfilled with the 
exception of the rate and torque limit during the very 
fast 3-segment slews. However, some operational 
obstacles like cat bed (thruster) heating [5] and 
communication were also adjusted in order to cope with 
this strategy. 



Settling Time 25% Margin 


Fig. 4 : Evolution of the angular momentum (L) in S/C frame along the Y-axis of the S/C, the angular velocity (w) and 
the resulting accumulated velocity increment ( AV ) along the required thrust direction over time (in minutes). 
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3.3 Global attitude 

As described below the AV WoL was performed around 

apocentre. During each orbit 7 off-loads of the angular 
momentum to +1-2 Nms along the S/C y-axis, and one 
extra off-load for each perilune (to compensate external 
torques) were performed. The perilune passage was 
commanded to be done Nadir pointing of the Z-axis to 
minimize the gravity gradient torque. Several Earth 
communication blocks where also inserted. The 
standard command generation software for the attitude 
and A L could be used. Only the interface with the 
optimisation software was adjusted to simulate a 
standard thruster behaviour. Starting on the 19/6/2006, 
overall more than 520 reaction wheel off-loads were 
executed during the first campaign. Later during a 
second campaign another 21 A L manoeuvres were 
executed. 

3.4 Performance 

The uncertainty in the off-load target was expected to be 
the largest error source due to the size of the off-load 
target threshold. However, over many off-loads these 
uncertainties averaged out as expected, resulting in only 
about 0.5-1% overall misperformance. The RCS 
modelling was supported by orbit determination and the 
respective performance calibrations. Further, the 
propagation of the reaction wheel was based on 
calibrated data from prior telemetry. The central 
problem was the estimate of the duration of an off-load. 
Since for performance reasons each slew needed to start 
as close as possible to the end of a A L manoeuvre, but 
still ensuring the finalization of the A L before this 
time, a 25% margin was assumed. A total AV of 1 1.8 
m/s using 2.5 kg Hydrazine to raise the pericentre by 90 
km in total was performed 

4. FULL OPTIMISATION SCHEME 

The objective of the optimisation was not to extend the 
mission or save up Hydrazine. The target was to reach 
height 0 for the first time in a point visible from the 
Earth and observable by the ground stations. 

Based on the small AV wol manoeuvres available, the 

optimisation strategy was to perform these in groups (of 
G, typically 5 or 7) distributed symmetrically with 
respect to the apolune (positioning them by true 

anomaly). Each such AV wol had the same magnitude 

and was pointing along the S/C velocity direction. For 
every orbit with manoeuvres, this would give an equal 
total impulse: 


AF orb =G-AF wol 


To gain enough degrees of freedom two of such series 
were planned one for end of June and one in July with a 
coast phase between them. A third slot in August was 
reserved for corrections. 

This full optimisation problem was solved numerically 
with MANTRA, an operational manoeuvre optimisation 
and trajectory propagation software tool. The most 
important perturbations considered are: 

• 3 rd Body, specially Earth and Sun 

• Moon W91 gravity field, up to 20x20 terms. 

• Solar radiation pressure and eclipse effect 

• Moon surface model W91 with a resolution of 
0.25 degrees in longitude and latitude. 

This approach had two main difficulties. It is heavily 
time consuming due to the high accuracy required for 
the long-term propagation. Further, the number of 
manoeuvres per orbit (G) and number of orbits per 
series (N) are discrete parameters that were input to the 
software package. These parameters must therefore be 
determined before hand manually: 

5. LINEAR APPROXIMATION 

The complete optimisation problem required an 
approximation, which allowed fast generation of the 
initial conditions, i.e. the numbers G and N. This was 
also of vital importance during operations to be able to 
generate the results on time. 


AV X AF. AV n 



Fig. 5 : Manoeuvre scheme 


5.1 Fundamental idea 

The low magnitude of the manoeuvres allowed 
assuming, that the effect of each is linear and 
independent of the previous one. So, a straightforward 
practical solution to the complete problem is to model it 
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as a superposition of single firings. The manoeuvres 
around the apolune are considered equivalent to a single 
instantaneous manoeuvre at the apolune [Fig. 5], In 
addition, the effect of the impulse at each apolune is 
independent of the previous one. Firing along the 
velocity direction the Gauss perturbation equations [8] 
show that only the in-plane parameters semi major axis 
(a) and eccentricity (e) -or equivalent the orbital period 
and perilune radius- will be affected. 

The change in these parameters is therefore modelled to 
be linear in the AV magnitude and for a succession of i 
apolunes (with manoeuvre) the approximation becomes: 


da 

dt 

de _ 
dt 


■r = 


_ 2-q 2 -V 

2 ■ (e + cos (l>)) 

V 

dco _ 2-sin(L>) 


2-a^Vl 


_y_ 

v. 


■7= 2-(e-l)~ 


dt 


e-V 


y = 0 


Here V is the velocity of the S/C, p the gravitational 
parameter of the moon, V the true anomaly and y the 
applied acceleration along the velocity. 


r pi ~ V + [ ' K i 

T,«T 0 +i-K 2 ' AV , 

where r pi is the perilune radius and T t the orbit period 

of the i-th orbit. K x and K 2 are proportional constants, 
related to the physics of the problem. 


The third equation shows that the argument of perilune 
does not change for a manoeuvre at apolune. The Earth- 
Sun perturbation keeps its value, and it is only changed 
by a small amount, resulting from the Earth and Moon 
movement. For the velocity in the apolune we have 



For N orbits, the approximation of the orbital period can 
be used to estimate the final time (taking into account 
that the first and last orbits count only half): 


J 1 N J 1 

t f = / 0 +— + 5 a - — - 

f 0 2 tT ' 2 

N 2 

= t n +N-T n + K, ■ AV 


where the following expressions where taken into 
account: 


r„ r„ n-k,- av 

2 2 2 


j]T r =N-T 0 +K 2 


f, A'-tA' + O 


■AV- 


Z' 


5.2 Perturbation equation 


Putting all together with the definition of the semilatus 
rectum, p = a ■ (l — e 2 ) we get for instantaneous 


increments 


Afl \-e_ AV_ 

a 1 + e V a 

AV 


Ae = -2 ■ (l - e ) 


V, 


by using the variables A and S i.e. 


r 1 + e 


£ = 


AV 


V„ 


the above equations become 


A a 

— = 2- A- 8 

a 

^ = 2-(l + X)-e 


In an intermediate step, we use the perturbation 
equations to prepare the deduction of a more rigorous 
estimate in the next section. 


As T is the osculating orbital period, then we get from 
the 2 nd Kepler law 


As the firing is tangential to the orbit, no out-of-plane 
parameter (inclination, line of nodes) is affected. 


AT 

T 


3 A a 


= 2- A-s 


It is known that the time derivatives of the semi major Therefore, the derivative of the radius of the perilune 

axis a, the eccentricity e and the argument of perilune fulfils 

CD for a tangential manoeuvre at the moon apocentre 
(U = 7T) are: 
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r p = r a -^1 

r a = const) 


Ar„ 


A/l 

~T 


2-(l + A)-s 


the error becomes greater or equal to £ 0 . In this case, 
the error becomes of the magnitude of AV orb . 


With these results, we can now determine the radius of 
the perilune and the period in each perilune. 


5.3 Constants and error estimation 

Here we give a quantitative estimate of the intuitive 
linearization idea from 5.1. 

First, the problem can be simplified by assuming that 
the increment in apolune velocity AV orb is small (of the 

order of 1 cm / S ) relative to the overall velocity at 
apolune V a (of the order of 1 km / S ), i.e. 


AV lcm/ s 

s 0 = * = 10 

V a0 1 km Is 


So the next approximations will keep only terms of the 
order of £ 0 . Second, the increase of £ due to the 
variation in the apolune velocity is small: 


V ai =v a0 +i- AV orb 


AV 


£, = 


orb 


> — ^ £ i — £ c\ ' 


AV, 


K 0 +i 'AV orb l + i-s 0 


orb 


a 0 



Fig. 6: Plot of A — A for several initial values of A 
and a value of £ 0 = 10 5 

With this the expressions for the perilune radius and the 


period after the i-th manoeuvre are 

r pt = r pn '(l + 2-(l + /lo)'s)' 

t, = r„.(i + 3. v«)' . 

Using the approximation, (l + x)” = 1 + n- x + o{x 2 
we can rewrite the previous equations as 


That limits the number of manoeuvres to be much 
smaller than l/<£ 0 =10 5 , and allows an easier 

comparison of terms with £ , since £ { and £ 0 are of the 

order of 10" 5 (of course depending on the magnitude 
of i) 

The next simplification is to keep the right hand term of 
the equations constant (see 5.2). To see how restrictive 
the number of orbits actually is, we study the evolution 
of the parameter /l using the full expression (but 
applying simplification of £ { ) and the simplified one: 

A = u + 2 ■ s ■ u • (i + u ) 

+ 2 • £ ■ A 0 • (l + A 0 ) . 


r Pi ar Po +2 ‘ / -( 1 + yl o )^ 

T i ~ Tq + 3 • i • Aq • £ 

These last expressions are again a linear approximation 
of the perilune and period as a function of the impulse. 
Moreover, a comparison with results from chapter 5.1 
shows: 


K, =2- , + A ° 


l 


0 


K 2 = 3 

V ' 2 v 

r a0 y aO 


So now, we see that for a series of less then 200 
manoeuvres our approximation has an error of the 

magnitude of A V orb . 


The SMART- 1 orbit had a high eccentricity with a 
typical value of 0.4 or higher. This results in a range for 

A of 0.3 < A < 0.5 . 

The graph [Fig. 6] shows the divergence between both 
series. The difference is lower than £ 0 while N<200 for 

the range of initial values considered, giving a tight 
limit for the number of manoeuvres. For bigger values, 


6. FURTHER CONSIDERATIONS 
6.1 Decay of the perilune 

With these results being valid for each manoeuvre series 
a model of the overall scheme is obtained. The plot of 
the osculating radius of the perilune at perilune as 
function of time allows us to obtain two powerful 
results. 
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Fig. 7: Perilune change per revolution up to SMART- 1 
impact (without manoeuvres). 


First and most important is that during the impact period 
the change in perilune radius must be strictly negative 
[Fig. 7], This will guarantee that the impact cannot be 
delayed by this effect, i.e. will happen in a near by orbit 
(max 5.5 hours error) 

The second result is that for long coast periods the 
decay of perilune can be approximated by its mean 
value (D) 


The equation for the final time is also more complex 
than the one of the final perilune radius, because the 
effect of the first will set the initial period for the second 
manoeuvre block. To pack the equation in a handy form 
we define the period at the end of each manoeuvre 

series as T a and T p respectively 

T a =T 0 +N a -K 2a - AFl 
T fi =T a +N fi -K 2/> - AFj 

^ a =N a -T 0 +^-K 2 a -AV a 

~ > 

N 2 

t f =t 0 + A T a +L a -T a + AT p + Lp -Tp . 

6.2 Other constraints 

Similar considerations for the illumination conditions 
and the visibility of the impact from Earth were made, 
to find initial values for the number of firing and free 
drift orbits N and L . 


We can combine the results (using the sub index a for 
the parameters of the first series, /? for the second 
series and L for the number of coast orbits) 


Ar =N -K -AV -D • L 


Arp — Np • Kp • AV p Dp ■ Lp 

r p f = V> + Ar a + Ar p 


For the evolution of the orbital period and the final time 
we must take into account that the Sun-Earth 
perturbation does not increment the energy of the orbit 
(as it is shown in the semi major axis plot [Fig. 8]) so 
the period is kept constant with a reasonable accuracy. 



Fig. 8: Evolution of the semi major axis up to SMART- 
1 impact (without manoeuvres). 


6.3 The final orbit 

From the equation of the final time, we cannot fully 
deduce the impact time since in case the impact is 
missed a complete new orbit must be added. A similar 
problem occurs if the impact happens one orbit before. 
There is also an error of about 1 minute between the 
true impact time and the perilune time since the impact 
happens in the descending arc of the orbit and not 
exactly at perilune. 

We can split the impact time variations in two parts of 
different magnitude: One big discrete correction of the 
order of an orbit period and a second one smaller and 
continuous of the order of the increment of period of a 
manoeuvre. For this second factor the following 
condition must be satisfied 

N a + L a + Np + Lp = const . 

7. OPERATIONAL RESULTS 
7.1 Iteration procedure 

The algorithm to find an optimal solution is as follows: 

a) Choose N a and N p 

b) Obtain L a and Lp with the linear model 
using a default value for the impulses 
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c) Implement a MANTRA case using these 
inputs. Set the optimisation constraints, 
including the height over the surface being 
greater than 0 for the last orbit and below 0 for 
the impact orbit. 

d) Run the case to obtain the manoeuvres 
magnitudes (A V a , A V p ) and the real impact 

time ^Impact = t{k = 0). 

e) Analyze the results for feasibility and stability, 
which includes studying the magnitudes of the 
manoeuvres to be close to valid values and 
ensuring an impact point visible from Earth. 

Once a first converged solution is found, one can 
transform it into initial conditions for further 
optimisations. This will result in shorter converging 
times. 

The first iteration was done to select the impact 
coordinates. The criteria to accept a case was only the 
visibility from Earth and the duration of the thrust phase 
for each orbit. For performance reasons, the starting 
point of the first A L must be below 120 degrees 
(measured in true anomaly from the apolune). 

With this data, the scientific community selected an 
impact time between 5:30 and 6:30 UTC on the 
3/9/2006 to be in the shadowed part of the Moon near 
the terminator, visible form from La Silla (Chile). 

The second iteration includes now the constraint for 
t impact • Several combinations of N a and N p were 
studied to select the final scheme. The final strategy 
chosen had 67 orbits each with 7 AV wol manoeuvres 

[Fig. 9/10]. It was split in 59 revolutions in June and 8 
in July. For August a correction slot was reserved to 
correct for possible misperformances. Since in average 
no such deviations occurred, it was not needed. 



Fig. 9: Evolution of the perilune radius (with 
manoeuvres) 



Fig. 10: Evolution of semi major axis (with 
manoeuvres) 


7.2 Impulse magnitude adjustments 

During operations, it is necessary to be able to adapt to 
the calibrated value of the manoeuvre performance. In 
order to quickly obtain a new approximation for the 

A V wol the following equations were used. 


' Pf 


' pf 


'pf 


= r p0 + N ■ M ■ K 2 ■ AV wol | 
= r p0 + N'-M'-K 2 . AC] 


N M AV 


= 


wol 


Pf 


N M AV 


= 1 


wol 


Starting with N and M and keeping the final perilune 
radius but changing the number of impulses or the 
number of orbits one can determine new values (N’ and 
M’). 


7.3 The Moon surface 

We used the model W91 for the lunar surface. This 
model is available and was obtained from the 
Clementine NASA mission [9]. It has a distribution of 
mean heights in a 0.25x0.25 degrees grid in latitude and 
longitude over the reference ellipsoid. The deviation of 
this model from the real Moon surface and the effect of 
the gravity potential for very low orbits were the main 
uncertainties. 


7.4 Last emergency manoeuvre 

On the 31st of August, 3 days prior to the impact, 
Antony C. Cook suggested a more accurate surface 
model. An analysis of the last orbit determination with 
this new model showed that the probability of hitting 
the edge of the Clausius crater one orbit before the 
predicted impact was too high [Fig. 11]. 
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Fig. 11: Study of the last orbit with the more accurate 
Antony C. Cook surface model [10]. 


To avoid the collision risk in an earlier orbit a correction 
manoeuvre was performed, to raise the perilune by 600 
meters. This ensured the S/C to be above the uncertainty 
in height of the A. C. Cook model. The nominal impact 



Fig. 12: Mosaic of images of the SMART- 1 impact site, 
obtained by the Advanced Moon Imaging Experiment 
(AMIE) on board SMART- 1, Mark R. Rosiek [11]. 


could be up linked. Thanks to the effort and experience 
of the operation team, the satellite could be recovered 
quickly and updated commands were up linked. 

8. CONCLUSION 

Using the on-board control logic, which converts a A L 
manoeuvre into a AV , and a flexible attitude generator 
the equivalent to a proper low thrust motor could be 
provided. The optimisation was based on the same 
software and similar ideas as for the earlier EP-firing 
phases. Clear interfaces in-between these two tasks were 
helpful to cut down the overall complexity. Although 
the feasibility was completely unknown at the time of 
the initial planning and many concerns arose, the final 
implementation proved to be very stable and smooth. 
Already with the first manoeuvre sequence sufficient 
accuracy was achieved, i.e. no further detailed 
refinement of the procedure was required, apart from 
standard calibrations. The science community was 
awaiting the impact on the morning of 3/9/2007. At 
5:42:21.759 UTC the SMART-1 carrier signal dropped 
to zero and the CFHT (Canada-France-Hawaii 
Telescope) infrared telescope recorded a flair at exactly 
the expected impact point (46.20W 43.44S) [Fig. 12]. 
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